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Abstract. We present the potential for discovering the Standard Model Higgs boson produced by the 



o : 

vector-boson fusion mechanism. We considered the decay of Higgs bosons into the W + W final state 



X 



with both VU-bosons subsequently decaying leptonically. The main background is tt with one or more jets 
produced. This study is based on a full simulation of the CMS detector, and up-to-date reconstruction 
codes. The result is that a signal of 5<r significance can be obtained with an integrated luminosity of 
12 — 72 fb^ 1 for Higgs boson masses between 130 < mn < 200 GeV. In addition, the major background 
can be measured directly to 7% from the data with an integrated luminosity of 30 In this study, we 

suggested a method to obtain information in Higgs mass using the transverse mass distributions. 

1 Introduction at 114.4 GeV for a 95% C.L. [T], and unitarity puts an 

upper limit of about 1 TeV. Even more constraining are 



One of the primary goals of CMS is to prove or disprove 
the existence of the Higgs boson. The LEP experiments 
set a lower limit on the Standard Model (SM) Higgs boson 
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the results of fits to precision electroweak measurements, 
which limit the mass of a Standard Model-like Higgs boson 
to be less than 194 GeV at 95% C.L.In extended Higgs 
sectors, there is often one scalar boson that resembles the 
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Fig. 1. Feynman diagram for Higgs boson production through 
Vector Boson Fusion. The Higgs boson decays into W's which 
further decay into electron / muon-neutrino pairs. 

Higgs boson of the Standard Model, and is responsible 
for electroweak symmetry-breaking. The mass of such a 
Higgs must also satisfy these constraints approximately. 
In the Minimal Supersymmetric extension of the Standard 
Model (MSSM), there is a more stringent bound coming 
from the internal constraints of the theory; the lightest 
Higgs boson must have a mass less than about 135 GeV. 
For these reasons, we focus on the mass region 120 < 
m H < 200 GeV. 

The two main decay modes of the Standard Model 
Higgs boson in this mass range are h — > bb and h — > 
W + W~. In the latter case, one of the W bosons may be 
off the mass shell. If the Higgs boson is heavier than about 
135 GeV, the WW* branching fraction will dominate, but 
it can be important for masses as low as 120 GeV. In 
this study, we consider the decay h — > WW* with the 
subsequent decay of the W^-bosons to two charged leptons. 
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Higgs bosons may be produced in pp collisions when ra- 
diated off the virtual VF-boson that is exchanged in the t- 
channel - this is called "Vector Boson Fusion" (VBF). The 
Feynman diagram for this process is shown in Fig. [1] This 
channel has good prospects for the discovery of a Stan- 
dard Model Higgs boson, especially if it is not too heavy 
because of the distinctive VBF topology which contains 
two jets with small angles with respect to the beam axis. 
Furthermore, when the Higgs decays to two VF-bosons, 
the presence of the hWW vertex both in production and 
decay of the Higgs boson gives a relatively clean determi- 
nation to the hWW coupling. Given the Higgs mass the 
Standard Model(SM) is completely determined, so that a 
measure of hWW coupling over-constrains the SM. This 
will be crucial to establishing the origin of electroweak 
symmetry breaking. 

The VBF mechanism was proposed as a potential dis- 
covery channel several years ago [3J. Our initial study of 
this channel for the CMS detector was carried out in 2002 [4] , 
with a number of simplifications. The conclusion of this 
previous CMS study was that a convincing signal for a 
Higgs boson with a mass of 120 GeV would be observed 
with about 70 fb _1 . In the present study, we repeat the en- 
tire analysis in the mass range 120-200 GeV, using the lat- 
est simulation and reconstruction software for CMS in or- 
der to verify and improve the 2002 study. A similar study 
of this channel for the ATLAS detector was performed in 
2004 using different generators and slightly different cuts 
0. 
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The VBF process is characterized by two forward jets 
with modest transverse momentum, Et ~ mw/2, sep- 
arated by a large rapidity difference. The Higgs boson 
signature is at low rapidity, with a pair of clean, iso- 
lated leptons and missing energy. The main backgrounds 
for this channel are the irreducible continuum W + W~ 
production, and tt in which both top quarks decay semi- 
leptonically. These backgrounds are particularly trouble- 
some when there are extra jets, j, in the event, so we have 
taken particular care with the generation of W + W~jj and 
ti events. 

2 Event Generation 

The signal process and the W + W~jj background have 
been simulated on the basis of a matrix-element calcu- 
lation using MadGraph [6]. For the tij background, we 
used the AlpGen [7] package which correctly simulates 
spin correlations. We simulated the parton showers us- 
ing Pythia [8]. MadGraph and AlpGen calculations are 
made leading order (LO). The parton distribution func- 
tions used by MadGraph and AlpGen are CTEQ6L1 and 
CTEQ5L1 respectively. The minimum transverse momen- 
tum cut on jets is 15 GeV, and the jet pseudo-rapidity 
is limited to < 5. We required a separation of any jet 
pair, namely, AR > 0.5, where AR = ^/(Ar]) 2 + (Acf)) 2 . 

Next-to-leading order (NLO) cross-sections differ from 
LO cross-sections by ~ 30% for a 120 GeV Higgs boson 
and ~ 10% for a 200 GeV Higgs boson [pj. However, since 
there are no NLO cross-section calculations for the back- 
grounds, the LO cross-sections are used consistently for 
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both signal and background processes in this study. The 
cross sections are listed in Table[TJ The 'electroweak' (EW) 
part of the W + W~jj process is defined as the subsample 
with no a s -dependent vertex in the diagrams, and the 
'QCD' part is the rest of this process. Note that the EW 
part is topologically very similar to the signal and hence 
is almost irreducible. 



Table 1. Production cross-section for the signal and main 
backgrounds 



Channel 


cross-section [pb] 


WW BR 


crxBR [pb] 


qqH m=120 


4.549 


0.133 


0.605 


qqH m=130 


4.060 


0.289 


1.173 


qqH m=140 


3.648 


0.486 


1.773 


qqH m=160 


3.011 


0.902 


2.715 


qqH m=180 


2.542 


0.935 


2.376 


qqH m=200 


2.177 


0.735 


1.600 


ttj 


736.5 


1. 


736.5 


WWjj QCD 


43.6 


1. 


43.6 


WWjj EW 


0.933 


1. 


0.933 



3 Detector Simulation and Event 
Reconstruction 

We processed the generated events through the CMS de- 
tector simulation software which is based on the Geant- 

4 simulation of the CMS detector. We simulated pile-up 
from out-of-time interactions representing the low lumi- 
nosity LHC running condition (~ 2 x 10 33 cm~ 2 s _1 ). Sub- 
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sequently, we processed digitized information (digis) using 
the CMS event reconstruction software. 

3.1 Trigger 

We refer to Ref. 10J for the presently planned trigger 
table. The inclusive single electron trigger has an Et- 
threshold of 26 GeV, which is too high for our purposes. 
Therefore we will augment this trigger with the di-electron 
trigger, which has a threshold of 12 GeV for both elec- 
trons. The pT-threshold for the inclusive single muon trig- 
ger is 19 GeV, which is well suited to this analysis. Con- 
cerning the e-/i channel, we plan to use the e+/i di-lcpton 
trigger, which will have a threshold of 10 GeV for each lep- 
ton. The efficiency for the Ll+HLT trigger with respect 
to our offline cuts varies from about 95% to 99% based on 
Ref. [11] . This presents no significant effect at the current 
state of our analysis. 

There will be lepton+jet triggers that should be very 
useful for this analysis if lower lcpton thresholds are needed. 
However, since the details for these triggers are not avail- 
able at this time, we have based our study solely on the 
leptonic triggers. 

3.2 Lepton Reconstruction and Identification 

We have used standard packages and selection criteria for 
muon and electron identification. Below, we describe our 
assessment of the identification efficiency. 
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3.2.1 Muons 

We use the "global" muon reconstruction, which takes 
muons found in the muon chambers (drift tubes, cath- 
ode strip chambers, and RPC's), and extrapolates them 
into the silicon tracker to pick up additional hits and bet- 
ter define the kinematics. This extrapolation takes into 
account the energy lost by the muon as well as multiple 
scattering. 

Muons are found within \rj\ < 2.4. The overall muon 
reconstruction efficiency in this angular range is ~ 95% 
for 10 < p T < 30 GeV and 97% for p T > 30 GeV. 

3.2.2 Electrons 

Electrons are reconstructed by combining super-clusters 
[T2Tfl~5] and Kalman tracks [H] • The track - super-cluster 
(SC) matching condition is AR < 0.15. Such tracks should 
have at least four hits, and transverse momentum pt > 
5 GeV. If several tracks satisfy these conditions, then the 
one having the least difference \pr — Et\ is taken. We 
reject the electron candidates if E T C < 10 GeV or \t]° \ > 
2.0. The probability for a generator level electron with 
Pt > 10 GeV and |r/| < 2.0 to be reconstructed within 
AR < 0.2 is - 92-98% for 10 < p T {gen) < 20 GeV and 
~ 98-99% for px > 20 GeV. These reconstructed electrons 
are said to be identified if they satisfy -Ehcal /-^ecal < 
0.05, |Z\?7(trk,SC)| < 0.005, E sc /p tlk > 0.8 and |l/£ sc - 
l/p uk \ < 0.06. 

An isolation variable is defined by taking the sum of 
the pt of all the tracks (except the electron candidate) 
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within a cone of AR SC < 0.2, and dividing by the E T C . 30% depending on the jet Et and r\ values. When applying 
The tracks entering this sum must have at least four hits, the correction to jets with \r\\ > 4, we used the correction 
Pt > 0.9 GeV, and |z trk — z e \ < 0.4 cm, where z is the parameters for the last interval \rj\ = 3.8 - 4.0. The polyno- 
position of the track along the beam line. We place the mial extrapolation is unreliable beyond pt = 200 GeV, so 
requirement that this isolation ratio be smaller than 0.2. we fixed the corrections above 200 GeV to those obtained 
The overall single electron efficiency for electron isolation at 200 GeV. The response to jets in the QCD di-jet sam- 
and identification is « 80% for 10 < pr < 30 GeV and pie is lower than the response to jets in the qqH sample. 
« 90% for pr > 30 GeV. The electron fake rate per jet This produces different correction functions. However, in 
is 3% for 10< p T <30 GeV and less than ss 0.1% for the current study, VBF tag jets are at high rj and have 
p T > 120 GeV calculated using the jets from W decay in at least pt > 30 GeV and for this part of phase space 
the associated production and using the forward jets in the differences between responses(or equivalently, the jet 
the qqH sample. correction functions) are very small. 

In the analysis, we used missing Et($t) calculated 
3.3 Jet and Missing E T Reconstruction and Correction from calorimeter hits. We corrected the^r using the sum 

of the Et difference between the corrected and uncor- 

The cell-level thresholds are set at least 2a above the noise 

rected jets for which the corrected jets have Et > 30 GeV. 

level to remove the effects of calorimeter noise fluctua- 
tions in jet reconstruction. This is important since we are 

mainly dealing with quite low-py jets in the current study. ^ Event Selection 

We reconstructed the jets using the "Iterative Cone" 

The strategy of the analysis is not complicated. We select 

algorithm, with a cone size of AR = 0.5 and a cone seed 

events with two forward jets separated by a large rapidity 

Et cut of 1 GeV. We removed the jets from an event 

difference, veto any event with additional central jets, and 

if they match a reconstructed electron within a cone of 

demand two energetic, isolated leptons in the central re- 

AR < 0.45. 

gion. Finally, we apply additional cuts on the kinematics 

We calibrated the reconstructed jets using the qqH sig- 

and the event topology. 

nal sample. Reconstructed jets arc first matched to gen- 
erator level jets within a cone of AR < 0.12. We fit the 

4.1 Forward Jet Tagging 

jet response to second-order polynomials as a function of 

generator-level jet Et for 20 different r\ regions covering The jets are ordered in Et after the corrections have been 
i] = to r] = 4 in bins of At] = 0.2. The difference between applied. The first two tag jets should be energetic, so we 
the corrected and uncorrected responses varies by 10% to require Eti > 50 GeV and Et2 > 30 GeV. Fig. [3] shows 
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Fig. 2. Arj = \rji — 772 1 distribution for the forward tagging jets which have Eti > 50 GeV and Et2 > 30 GeV for a) qqH, 
rriH = 120 GeV and backgrounds b)tij, c) EW WWjj and d) QCD WWjj. Note that the EW WWjj background is basically 
irreducible. 

the rapidity separation \Arj\ between these two most en- 4.2 Central Jet Veto 

ergetic jets, for the signal(a) and the backgrounds (b-d). It 

is clear that the jets for signal events are well separated in In the signal process, there is no color exchange between 
rapidity, and we apply the cut \Ar)\ > 4.2. We also make the protons, and consequently any additional jets will tend 
sure that they fall in opposite laboratory hemispheres by to be radiated in the forward direction. In contrast, the 
requiring 771 • J72 < 0. backgrounds will tend to have additional jets in the cen- 

tral region, especially the tij process. We take advantage 
of this distinction by vetoing events with additional jets in 
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Fig. 3. r/o = 773 — (771 + 772 ) /2 for the third jet. 77 of the third jet with respect to the average of the two forward jets. For signal 
a) qqH, m H = 120 GeV and backgrounds b) tij, c) EW WWjj and d) QCD WWjj. 



the central region. In particular, we consider any jet with 
Et3 > 20 GeV and compute the rapidity with respect to 
the average of the two forward jets: 770 = 7/3 — (771 + 772 )/2. 
We veto the event if \rjo\ < 2. See Fig. [3] for distributions 
of both signal and background. The probability to find a 
fake jet from pile-up events for low luminosity LHC run- 
ning is shown in Fig. 0] as a function of the Ex threshold 
for the central jet veto. The fake rate is defined as the rate 



for pile-up jets satisfying the central jet veto condition in 
an event where there are no real jets satisfying those con- 
ditions. Therefore, the fake rate is just the rate of events 
mistakenly rejected due to pile-up. The loss of events for 
a E T threshold of 20 GeV is only about 2%. 

The effect of the Et threshold for the central jets on 
the final cross sections and significances for the 120 GeV 
signal and for the background are displayed in Fig.0 Here, 
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Fig. 5. The effect of the Et threshold for the central jet veto. For a) signal evernts, qqH with mn = 120 GeV and background 
events b) ttj. c) the S/B ratio and d) the significance for a 30 /6 _1 integrated luminosity. 



the significance is defined as S/wB, where 5* and B rep- lepton must be low since one of the two W's in the Higgs 
resent the numbers of signal and background events. decay is off the mass shell for low Higgs masses. Fig. |6] 

shows the px spectra for electrons in the signal process 
4.3 Lepton Kinematics (Mjj = 120 GeV). We reject events with more than two 

leptons. The two leptons must be well separated from all 



We require two opposite-sign leptons in an event. The 
most energetic lepton must have pri > 20 GeV, and the 
other, pT2 > 10 GeV. The pr-threshold for the second 



jets with ARij > 0.7. 
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On > 26 GeV AND p T2 > 10 GeV) 
OR 

(pti > 20 GeV AND p T2 > 12 GeV) . 
Since the leptons come from the W's that come from 

the centrally-produced Higgs boson, we require them to be 

central. If rj^i is the forward-tag jet having higher-rapidity, 

and rji is that of the lower-rapidity forward-tag jet, then 

our requirement can be written t]\ + 0.6 < r\i < ?7hi — 0.6. 

This condition must be satisfied by both leptons. Fig. [7] 

Fie. 4. Fake central jets fraction per event as a function of , , •, r , , i , j ; / 

snows the distributions of the related quantity, rj e = (rje — 
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at least one jet(due to pile-up) satisfying the central jet veto 
conditions, with no "real" jets satisfying the central jet veto 
condition in that event. 



+ ? ?2)/2) x 4.2/Arj. This quantity is sensitive to the 
•q distribution of leptons with respect to the forward tag 
jets. 




Higher P T Electron 
Lower P T Electron 



200 250 

P T (GeV) 



Fig. 6. Electron pr spectra for the signal process, qqH, when 
m H = 120 GeV 



In light of the thresholds for the electron triggers, we 



4.4 Further Kinematic Requirements 

After the forward-jet tag, the central jet veto, and the 
lepton kinematics cuts, we are left with a sample which 
still has a large contamination from background processes. 
We can further reduce this contamination with some ad- 
ditional kinematic cuts. 

First, we require the di-jet mass to be greater than 600 GeV 
(see Fig- [HI) ■ Next, we look at the overall pr-balance in the 
event, by computing the vector sum of the transverse mo- 
menta of the two leading jets, the leptons, and the missing 
energy. The magnitude of that sum should be less than 
40 GeV (see Fig.©. 

When it comes to the leptons, we require a di-lepton 



modified our pr requirements slightly in the di-clcctron mass M« < 80 GeV (see Fig. [TU)) . This value is lower 
channel. An event is selected if it has two electrons which than the Z-mass, so that leptonic Z-decays do not affect 



satisfy: 



the current analysis. A useful distinction arises in the rel- 
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Fig. 7. Centrality of the leptons, using the quantity r\ t defined in the text for a) qqH, tuh = 120 GeV and backgrounds b) ttj, 
c) EW WWjj and d) QCD WWjj. 



ative azimuthal angle of the two leptons due to the spin-0 
nature of the Higgs boson (see Fig. [TTJ [T^|) . We take ad- 
vantage of this discriminant and require A<j) < 2.4 radians. 
Finally, we require that the "WW transverse-mass" be not 
too high when looking for Higgs bosons with mass below 
150 GeV. The cut is that M T ,ww < 125 GeV, where 
M TWW = yJ(ff T +PT,a) 2 - Wr + Ptm) 2 - See Fig. [TJ 
and [14] for distributions of this quantity. 



4.5 Additional Cuts 

Additional cuts may be required in order that bbjj and 
rrjj backgrounds not pose a problem. The additional cuts 
57.3Z\</>(«,$ T ) + 1.5p™ sgs > 180 and 12 x 57.29A£(«, / 
E T ) +Pt SSS > 360 (where A(j){U,E T ) is in radians and 
p T isss is in GeV units), and also^V > 30 GeV if p™ sgs < 
50 GeV, are imported from Ref. Here, py lggs is the 
vector sum of the transverse energy of tag jets. The dis- 
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Fig. 8. Invariant mass distributions for the two forward tag jets, for a) qqH, ttih = 120 GeV and backgrounds b) ttj, c) EW 
WW}} and d) QCD WW}}. 



tribution of signal events in the A(/)(££,I$t)-Pt 1SSS plane is 
displayed in Fig. fT51 



The Drell-Yan production of di-lepton pairs, 7* — > 
£ + £~, has a large cross-section. In order to reduce this 
background sufficiently, we impose a di-lepton mass cut 
Mu > 10 GeV and we require ty T > 30 GeV when the 
leptons have the same flavor (see Ref. [3]). 



Finally, we impose the cut A(j>{££,$ T ) + A<p(££) < 3 
radians, which increases the signal-to-background ratio. 
Fig. [16] shows distributions of this quantity. The resolu- 
tion of the quantity Acp(££,fix) 1S improved by the $t 
correction. The additional cuts imposed after the trans- 
verse mass cut were determined for generator level anal- 
ysis. Therefore, we did not include these cuts in the sig- 
nificance, background or mass estimation and their effect 
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Fig. 9. The overall pr-balance in the event. See the text for 
tij, c) EW WWjj and d) QCD WW}}. 

is seperately shown in Tables 4-6. Work is in progress to 
confirm their effect after full detector simulations. 

5 Results 

The total accepted signal cross-sections range from about 
0.8 fb up to 7.2 fb, depending on the Higgs mass. They are 
listed in Table O The contributions from the e + e~ and 
/i + n~ channel are very similar, and the e fi^ channels 




Pt,o,a,[GeV] 



explanation, for a) qqH, ran = 120 GeV and backgrounds b) 

are twice as large due to branching ratios. The total effi- 
ciency is 3-6%, depending on run- The background cross- 
sections are somewhat larger, and there are two back- 
ground values corresponding to the "low-mass" and the 
"high-mass" cuts - see Table [2j 

We computed the significance S c p of an excess of events 
over the tij and W + W~ jj backgrounds, assuming an in- 
tegrated luminosity of C — 10, 30 and 100 fb -1 . S c p 
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Fig. 10. Di-lepton invariant mass distribution after jet and lepton cuts, for a) qqH, mn = 120 GeV and backgrounds bttj, c) 
EW WWjj and d) QCD WWjj. 



is the probability calculated assuming a Poisson distri- 
bution with Nb background events to observe equal or 
greater than a total number of signal and background 
events {N$ + Nb), converted to an equivalent number of 
sigmas for a Gaussian distribution [15] . The code to cal- 
culate S c p is taken from Ref. [IB] , 

The background uncertainty is included in the calcula- 
tion. This uncertainty comes from the statistical error in 



the background estimation and amounts to about 12% at 
10 ftT 1 , 7% at 30 ftT 1 and 4% at 100 fb" 1 . See SectionO 
for a discussion of the background estimation. 

The results are summarized in Table [3J Even for a 
Higgs mass as low as 130 GeV, a 5er signal can be ob- 
tained with a reasonable integrated luminosity. For higher 
Higgs masses, a very strong signal would be expected, 
and prospects for a measurement of the cross section for 
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Fig. 11. The distribution of the difference in azimuthal angle between the two leptons, A<fi after jet and lepton cuts, for a) signal 
events, qgH, ran — 120 GeV and backgrounds bttj, c) EW WWjj and d) QCD WW]]. 



pp — » qqH become more promising. Fig. 1 171 shows the sig- 
nificance for an integrated luminosity of 30 fb _1 as a func- 
tion of run, and Fig. ll8l shows the minimum integrated lu- 
minosity needed for a 5a signal also as a function of ran. 
The individual cut efficiencies with respect to the starting 
cross-section for 120 and 160 GeV Higgs bosons and the 
backgrounds are shown in Tables I4I5I6I for each channel. 



Concerning systematics, we have first considered the 
impact of the jet energy scale. The expected jet energy 
scale uncertainty in CMS is about 3%. For the ttj back- 
ground the scale uncertainty after correction is about 5% 
for px > 30 GeV. In this analysis, the two tag jets are re- 
quired to have EpTi >50 GeV and Eti >30 GeV and we 
reject additional jets in the central region if their Et > 20 
GeV. For the jets with Et ~ 20 GeV, the cross-section un- 
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Table 4. Accepted signal (for niH = 120,160 GeV) and major background cross-sections in fb for the eevv final state. 



Cut 


qqH 


qqH 


ttj 


WWjj 


WWjj 




120 


160 




EW 


QCD 




5.261 


26.97 


8617. 


10.74 


514.3 


Eti > 50,E T2 > 30 GeV 


3.742 


18.70 


6743. 


8.838 


296.4 


An > 4.2 


1.217 


6.067 


184.2 


2.195 


12.22 


r?i x r)2 < 


1.215 


6.054 


183.1 


2.193 


12.18 


Mjj > 600 GeV 


1.073 


5.367 


147.2 


2.071 


9.052 


Pt — balance cut 


0.653 


3.353 


54.89 


1.021 


3.298 


Central Jet Veto 


0.401 


2.309 


15.04 


0.631 


1.490 


> 2 good leptons w opp. charge 


0.269 


1.915 


10.98 


0.483 


0.695 


p T > 20, 10 or p T > 26, 12 GeV 


0.250 


1.838 


10.59 


0.475 


0.675 


\AR(j,l)\ >0.7 


0.250 


1.830 


10.33 


0.471 


0.662 


Req. leptons between jets 


0.235 


1.712 


4.990 


0.417 


0.430 


Mu < 80 GeV 


0.235 


1.683 


2.386 


0.144 


0.205 


A<t>u < 2.4 


0.220 


1.587 


2.088 


0.127 


0.192 


M T ,ww < 125 GeV 


0.183 




1.139 


0.081 


0.093 


57. 29 A(j){ll,B T ) + 1.5p T (H) > 180 & 












12 x 57. 29 A(j>(ll,^ T ) +Pt{H) > 360 


0.161 




0.936 


0.069 


0.073 


M ;i > IO&^t > 30(ee, aa) 


0.115 




0.800 


0.053 


0.060 


A<t>(ll,B- T ) + A<t> a < 3 


0.090 




0.420 


0.031 


0.033 


High Mass Cuts 












No M T ,ww Cut 




1.587 


2.088 


0.127 


0.192 


57. 29 A4>(ll,tf T ) + l-5p T (H) > 180 & 












12 x 57. 29 A(j>{ll,^ T ) +Pt{H) > 360 




1.501 


1.885 


0.114 


0.172 


Mu > 10 GeV&B'T > 30(ee, w ) GeV 




1.303 


1.736 


0.098 


0.152 


A<j>(ll$ T ) + Acj>u < 3 




0.862 


0.651 


0.052 


0.046 
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Table 5. Accepted signal (for niff = 120,160 GeV) and major background cross-sections in fb for the e\ivv final state. 



Cut 


qqH 


qqH 


tti 


WWii 

V V V V JJ 


WWii 




120 


160 




EW 


QCD 




1U.0 r 


0o.z4 


1 TO OA 


Ol AO 

zl.4o 


i non 
lUzy. 


it 1 t? ---- on ~\ j 
h/Tl > 0U,£/T2 > oU Lrev 


t oon 
/ .zyU 


O0.04 


loozU. 


1 *7 OO 
1 / .ZZ 


OO 1 .1 


LsTj .> 4.Z 


Z.40o 


1 9 

1Z.00 


OOo.O 


4.0OO 


9/1 *3G 

z4.oy 


T)l X Tj2 < 


O A C /I 

z.4o4 


i o c c 
12.00 


orr r 

ooo.o 


A C OiS 
4.0ZD 


O A OK 

z4.zo 


Mjj > oUU Lrev 


o i /in 
z\14y 


1 1 no 
11. Uo 


zoz.U 


/i oon 


1 Q OQ 


Pt — balance cut 


l.oyo 


/ .oylt 


117/1 
1 1 ( .4 


o /I nc; 
Z.4U0 


fi OQ7 

o.zo r 


central Jet Veto 


U.o < y 


C 1 OQ 

O.lzo 


QO 7H 

oz. / U 


i cno 
l.OUz 


/I 1 OQ 

4.1zo 


> 2 good leptons w opp. charge 


U.O 


A OOO 
4.000 


O C H7 


1 1 Q£! 

l.loO 


o i no 
Z.lUz 


Pt > 20, 10 GeV 


0.544 


4.079 


23.47 


1.131 


1.975 




u.ooy 


A 0^9 


91 71 
Zl. / 1 


i i nn 


1 8R1 


R,eq. leptons between jets 


U.OUD 


^ 7/1 & 


i n fin 


u.yzu 


i nfiR 

1 .UDo 


iviii < oU Ljev 


U.OUO 


o.uoO 


0.U14 


U.oUl 


U.44 ( 


A(f)n < 2.4 


0.480 


3.497 


4.216 


0.245 


0.394 


MT^W \ ±ZO V 


n /inn 




9 fi91 
Z.UZ 1 


n i aa 

U. l^i^i 


0. 907 
U.ZU ( 


nn /i i ( n ~r~f \ i -i r* / rr\ - -t or\ 

57.29A4>(U,Pt) + 1.5pt(H) > 180 & 














0.329 




1.880 


0.109 


0.153 


i^T ^ OU vie V 11 pT\±l ) OU Vjrc V 


u.ozo 




1 89^ 


U. 1UO 


U. lOO 


AMU fid 4- A^,, s ^ 


u.zoy 




n 7QS 


U.UDO 


U.Uo 


High Mass Cuts 












No M T ,ww Cut 




3.497 


4.216 


0.245 


0.394 


57. 29 A<j>(ll,$ T ) + l-5p T (H) > 180 & 












12 x 57. 29 A(f>(ll,tf T ) + Pt{H) > 360 




3.105 


3.418 


0.202 


0.334 


#r > 30 GeV if p T (H) < 50 GeV 




3.084 


3.361 


0.199 


0.334 


A4>{ll,tf T ) + A<j> n < 3 




2.003 


1.709 


0.107 


0.173 



E. Yazgan, J. Damgov, et al.: Search for a SM Higgs Boson in CMS via VBF in the H^WW- > Ivlv Channel 
Table 6. Accepted signal (for mu —120,160 GeV) and major background cross-sections in fb for the \x\xvv final state. 



Cut 


qqH 


qqH 


tti 


WWii 


WWii 




120 


160 




EW 


QCD 




o.loo 


on /i /i 
zy.44 


ool (. 


11). / / 


K 1 O T 

Olz. / 


it 1 x c v on r i „"\ j 
£jt\ > 0U,ilT2 > oU Oev 


Q Q CT 




OOzl. 


C QQO 

o.ooz 


OQO CC 

zoz.O 


LsTj .> 4.Z 


1 971 
1 .Z ( 1 


7 QQ1 

( .oyi 


1 7S n 

1 ( o.U 


9 

Z.oDO 


1Z.11 


T)l X Tj2 < 


1 0£! O 
l.ZOO 


T QTK 

(.3(0 


1 TC? T 

1 ro. / 






Mjj > DUU Kjev 


i i no 


a coo 
O.Ozz 


1 QA 7 

toy. / 


O OC1 

z.zOl 


C ACQ 

o.yoo 


Pt — balance cut 


U.o04 


/I G/17 


00. / 


l.OoO 


K 7^Q 
O. / DO 


central Jet Veto 


U.OOz 


Q COQ 
O.OZO 


ly.Oo 


l.UU / 


o.loy 


> 2 good leptons w opp. charge 


n A QA 


O Qfll 

z.oyl 


1 1 1 

lu. 11 


n 770 

U. / / z 


1 A TO 

1.4 f z 


Pt > 20, 10 GeV 


0.327 


2.605 


14.30 


0.716 


1.324 


1 A T?( i JM \ 07 
l^-^Ui L )\ > u -' 


u.oiy 


9 ^7 


1 1 

ii. oy 


U.uoU 


i im 

l.loD 


R,eq. leptons between jets 


u.zyu 


9 9QS 
z.zyo 


0.4D1 


U.OOD 


U.04O 


jv^n ou ^re v 


n 9GH 

u.zyu 


9 99fi 
Z.ZZO 


9 Q71 
Z.O ( 1 


n i on 
u.iyu 


n 971 

U.Z ( 1 


A(j)u < 2.4 


0.273 


2.102 


2.024 


0.165 


0.252 


MT^W \ ±ZD V 


U.ZOo 




l.UDO 


n no9 
u.uyz 


n 1 1 q 


nn A i ( n ~r~f \ i -if / rr\ - -t or\ D 

57.29A4>(11,Pt) + 1.5pt(H) > 180 & 














0.200 




0.826 


0.075 


0.095 










u.uuu 


n 07^ 




A 

U.±o4 




U.4ZO 


U.UOl 


n n^9 

U.UDZ 


High Mass Cuts 












No M T ,ww Cut 




2.102 


2.024 


0.165 


0.252 


57. 29 A<j>(ll,$ T ) + 1.5p T {H) > 180 & 












12 x 57. 29 A(j>{ll,tf T ) +Pt{H) > 360 




1.908 


1.785 


0.147 


0.229 


> 30 GeV if p T (H) < 50 GeV 




1.681 


1.678 


0.132 


0.205 


A4>{ll,tf T ) + A<j> n < 3 




1.229 


0.746 


0.092 


0.119 
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Fig. 13. The transverse mass of the two W bosons, Mt,ww , for a) signal events, qqH, inn = 120 GeV and backgrounds b) ttj, 
c) EW WWjj and d) QCD WWjj. 

certainty after jet correction is about 10%. We re-computed 10% increases the significance by 0.3 - 3.4% depending 
all yields after scaling the raw jet energies up and down on ran ■ This is a systematic uncertainty on the signal 
by 10%. In general, signal and background yields correlate, cross section, 
so the impact on the significance with a 10% jet energy 

scale uncertainty is less than - 8 - 10% at 30 We also used the p y thia event generator for our signal 

as an alternative to MadGraph. For mn — 120 GeV, the 
We also tested our results for the significances to errors significance obtained with Pythia is higher by 30% for a 
in thei^T scale. Increasing the^r scale by 10% decreases luminosity of 100 fb~ , while for mn = 160 GeV, it is 
the significance by 9 - 11%. Decreasing the $t scale by higher by 10 %. 
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■VwtGeV) Mr,ww(GeV) 

Fig. 14. The transverse mass, Mt,ww , distributions for signal and background, with Higgs mass = 120, 130, HO, 160, 180 
and 200 GeV respectively shown in a),b),c),d),e),f). The Lower plot (light grey) is the signal, the middle plot(dark grey) is the 



background, and the black histogram is the sum. 
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Fig. 16. a) A(f>(U,I$T) vs. A<f>u and b) the som of Acf>u and A(f>(ll,I^T) for qqH tuh = 120 and for ttj bacground c) and d). The 
lines correspond to A<f>u + A<\>(\,IJ$t) = 3 radians. 



We found that the production cross-section depends 
on the choice of scale (renormalization scale x factorization 
scale) for the ttj background. The ttj cross-section is 736.5 
pb as reported in Table 1, with the definition of the scale 
imy, where m\ = rntop + Pt an d the sum is over final 
state light partons. However, if we change the definition 
of the above sum to include all the final state partons in- 
cluding the heavy quarks, then the cross-section decreases 
to 530 pb. These two definitions of scale are the defaults 
in AlpGen 1.3.3 and 2.0.x respectively. We found that the 
choice of scale does not affect the kinematics of ttj at 



all. Moreover, the cross-section and kinematics of the qqH 
process are not affected by the choice of scale. The signif- 
icance with the new scale choice is ~ 18% higher. There- 
fore, the uncertainties in the computed ttj background 
make it very important to measure the background di- 
rectly in the experiment. 

It should be pointed out that the statistical significance 
of our analysis is generally a factor of ~ 2.6-3.2 lower 
than the significance reported in the study for the ATLAS 
detector [5]. There are several reasons for this difference. 
First of all, the ttj cross-section used in Ref [S] is smaller 
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Fig. 12. The A<f> distribution between the two leptons after jet 
and lepton cuts for qqH, mu = 200 Ge V 




Table 2. Summary of accepted cross sections, in fb. A se- 
ries of assumed Higgs boson masses is shown, as well as the 
backgrounds for the "low-mass" and "high-mass" cuts. 

accepted cross-sections (fb) 
channel e+e~ e+fi+ fi+ jj,~ sum 



3 

A«ll#r) 



Fig. 15. The azimuthal angle difference in radians between the 
dilepton momentum vs. the missing Et vs p T T * 9gs for qqH with 
run = 120 GeV. The lines correspond to the cuts 57.29A<f>(ll, / 
E T ) + 1.5p^ iBga > 180, 12 x 57.29A£(H,-9r) + P?" s > 360. 

than the cross-section we use by about a factor of 0.7. 
Furthermore, the ATLAS study includes the gluon-gluon 
fusion channel for Higgs production which increases the 
signal by about 10%. 



"low" mass 

qqH, m H = 120 GeV 0.183 0.400 0.253 0.836 

qqH, m H = 130 GeV 0.387 0.854 0.601 1.842 

qqH, m H = 140 GeV 0.617 1.341 0.955 2.913 

ttj 1.139 2.621 1.065 4.825 

W+W-jj (EWK) 0.081 0.144 0.092 0.317 

W + W~jj (QCD) 0.093 0.207 0.119 0.419 

all backgrounds 5.561 





"high" mass 






qqH, m H = 160 GeV 


1.587 


3.497 


2.102 


7.186 


qqH, m H = 180 GeV 


1.362 


3.089 


1.837 


6.288 


qqH, m H = 200 GeV 


0.815 


1.703 


1.087 


3.605 


tij 


2.088 


4.216 


2.024 


8.328 


W+W-jj (EWK) 


0.127 


0.245 


0.165 


0.537 


W+W-jj (QCD) 


0.192 


0.394 


0.252 


0.838 


all backgrounds 








9.703 



Another important difference between the two analy- 
ses concerns the central jet veto. Our signal simulation 
generates a larger number of central jets compared to 
the ATLAS study, which used the PYTHIA Monte Carlo 
event generator. When we compare the signal efficiency 
after all cuts using PYTHIA instead of MadGraph, we 
find a difference of ~ 5 — 50%. Finally, the very definition 
of significance (S c p) differs between the two studies. The 
ATLAS study used a definition which gives a value which 
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Fig. 17. Significance of the Higgs signal as a function of Higgs 
mass for a 30 fb~ 1 integrated luminosity. 



Table 3. Significance of an excess as a function of Higgs 
mass, for three assumed integrated luminosities. The last col- 
umn shows the minimum luminosity required for a 5<r excess. 

Higgs mass significance £-min 

(GeV) 10 ftT 1 30 fbT 1 100 fbr 1 (fb _1 ) 

120 
130 
140 
160 
180 
200 



0.72 


1.35 


2.60 


340 


1.77 


3.04 


5.85 


72 


2.68 


4.79 


8.33 


33 


4.54 


7.00 


13.0 


12 


3.95 


6.22 


11.6 


15 


2.31 


4.03 


6.99 


45 
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Fig. 18. Minimum integrated luminosity (fb -1 ) needed to 
obtain a 5a excess over the ttj + W + W~jj background as a 
function of the Higgs mass. 

is ~ 9-14% higher for the same number of signal and 
background events. If the number of background events 
is reduced, the apparent improvement in the significance 
increases more dramatically than for our measure of signif- 
icance. Thus the uncertainty of ~ 9-14% should be taken 
as a lower limit for this particular factor. Considering all 
of the above, the differences between our results and those 
reported in Ref. [5] can be understood. Nonetheless, these 



considerations show that there still are uncertainties in the 
modelling of this channel which should be investigated by 
both experiments. 

5.1 Background Estimation from the Data 

For the Higgs masses considered here, there is practically 
no signal with M u > 110 GeV - see Fig. [101 For the 
present discussion wc define this as the signal-free region. 
Fig. [TTJ] shows the Mu distribution computed with looser 
cuts (no central jet veto, no pT-balancing cut, \ Arj\ > 3.5, 
T^io + 0.3 < rjg < ?/hi — 0.3) and the full analysis cuts. 
The number of events with M« > 110 GeV is desig- 
nated by "a" for the distribution with looser cuts and 
by "c" for the full analysis cuts. The number of events for 
Mu < 80 GeV is designated by "6" for the distribution 
with looser cuts and by "d" for the full analysis cuts. The 
region 80 < Mu < 110 GeV is excluded from the cal- 
culation in order to avoid any background coming from 
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Fig. 19. The Mu distribution computed with looser cuts and 
full analysis cuts. 

Z — ► £ + £~. Since Mu > 110 GeV represents the signal- 
free region, we can use the numbers a, c and b to estimate 
the number of background events in the region where we 
expect the signal (i.e., d). Using the simulations, we find 
that c/a = 0.097 and d/b = 0.098. The error on this esti- 
mation is dominated by the statistical uncertainty which 
is \/c/c w 7%. In order to obtain the background distri- 
bution in Mt^wWj we take the distribution obtained with 
the looser cuts and scale it by the factor of 0.098. A com- 
parison of the real and rescaled background distributions 
is given in Fig. [20l which indicates that this "data driven" 
works quite well. 

5.2 Sensitivity to the Higgs Mass 

The above significance estimates are for a pure "counting 
experiment" . We can, in addition, use the information con- 
tained in the distribution of Mt,ww with regard to the 
Higgs mass. We infer the mass of the Higgs boson from the 
observed distribution in Mt,ww by subtracting the data- 




Fig. 20. The transverse mass, Mt,ww, distribution for esti- 
mated(dashed) and real(solid) background. 

driven estimate of the background Mt,ww distribution 
from the distribution obtained with the full set of analysis 
cuts. The estimated and real Mxyvw distributions for sig- 
nal events are shown in Fig. [21] for several different Higgs 
boson masses. The inferred and the real mean values and 
shapes approximately agree. 

In an effort to obtain a quantitative measure of m#, 
we can use signal Mf.ww distributions as templates to be 
compared to the observed distribution. The comparison is 
done using the Kolmogorov-Smirnov test, and the results 
are shown in Fig. [5U A value close to one indicates a 
good match between the shapes. Comparing the means 
and shapes of the observed and template distributions, 
we can differentiate between Higgs boson masses for the 
cases of 160, 180 and 200 GeV, and for low masses (120 
- 140 GeV). To differentiate between the cases of 120, 
130 and 140 GeV Higgs mass, we must reduce the tij 
background more or we must have data corresponding to 
an integrated luminosity greater than 50 fb _1 . 
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Fig. 21. Estimated(dashed) and real(solid) Mt,ww distributions for signal events, with Higgs mass of 120,130,140,160,180 and 
200 GeV shown in a),b),c),d),e) and f) respectively. 
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Fig. 22. Kolmogorov test function for estimating the Higgs boson mass for Higgs masses of 120,130,14-0,160,180 and 200 GeV 
shown in a),b),c),d),e) and f) respectively. 

6 Conclusions boson fusion channel. We utilize the final state in which 

both W bosons decay to electrons or muons. Our study 
We have presented an analysis meant to isolate a discov- is based on a ful1 simulation of the CMS detector and an 
ery signal for a Standard Model Higgs boson in the vector- up-to-date version of the reconstruction codes. Further- 
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more, we have generated the main backgrounds, ttj and 5. S. Asai et al. Prospects for the search for a standard model 

W + W~jj, as accurately as is presently possible. Higgs boson in ATLAS using vector boson fusion, Eur. Phys. 

The results of our study are encouraging, and indi- J- C32S2 (2004) 19 

cate that an excess signal with a statistical significance 6 - F - Maltoni and T. Stelzer, MadEvent: automatic event gen- 

of over 5a can be obtained with an integrated luminos- eration with MadGraph, JHEP 0302 (2003) 027 

., r -,-,„_! , „ n „ _i „ TT . ,, 7. M.L. Mangano, M. Moretti, F. Piccinini, R. Pittau, A. 

lty of > 11 fb and < 72 fb tor Higgs masses m the 

Polosa, ALPGEN, a generator for hard multiparton processes 

range 130 < run < 200 GeV. Our analysis also shows that 
the background can be measured to 7% accuracy directly 
from the data. This uncertainty is dominated by statistics 
for 30 fb -1 . Finally, we suggest a method to obtain infor- 
mation on the Higgs mass using the shape of the Mt,ww 
distribution. 



in hadronic collisions, JHEP 0307 (2003) 001 
8. Torbjorn Sjostrand, Leif Lonnblad, Stephen Mrenna, 



PYTHIA 6.2: Physics and Manual, hep-ph/0108264 
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